Bertrand spacetimes as Kepler /oscillator potentials 



Angel Ballesteros a , Alberto Enciso b , Francisco J. Herranz c and Orlando Ragnisco d 

a Depto. de Fisica, Facultad de Ciencias, Universidad de Burgos, 09001 Burgos, Spain 
E-mail: angelb@ubu.es 

b Depto. de Fisica Teorica II, Universidad Complutense, 28040 Madrid, Spain 
E-mail: aenciso@fis.ucm.es 

c Depto. de Fisica, Escuela Politecnica Superior, Universidad de Burgos, 09001 Burgos, 
Spain 

E-mail: fjherranz@ubu.es 

d Dipartimento di Fisica, Universita di Roma Tre and Instituto Nazionale di Fisica 
Nucleare sezione di Roma Tre, Via Vasca Navale 84, 00146 Roma, Italy 
E-mail: ragnisco@fis.uniroma3.it 



Abstract 

Perlick's classification of (3 + l)-dimensional spherically symmetric and static 
spacetimes ( A4 , r] = — y di 2 + g) for which the classical Bertrand theorem holds 
[Perlick V 1992 Class. Quantum Grav. 9 1009] is revisited. For any Bertrand 
spacetime (A4,rj) the term V(r) is proven to be either the intrinsic Kepler- 
Coulomb or the harmonic oscillator potential on its associated Riemannian 3- 
manifold (M,g). Among the latter 3-spaces (M,g) we explicitly identify the 
three classical Riemannian spaces of constant curvature, a generalization of a 
Darboux space and the Iwai-Katayama spaces generalizing the MIC-Kepler and 
Taub-NUT problems. The key dynamical role played by the Kepler and oscillator 
potentials in Euclidean space is thus extended to a wide class of 3-dimensional 
curved spaces. 
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1 Introduction 



It is a well known, but nonetheless quite remarkable fact that all the bounded trajec- 
tories in Euclidean space associated with the Kepler-Coulomb (KC) potential V(r) = 
a/r + (3 are periodic, so that they yield close orbits, i.e. ellipses. This is in striking con- 
trast e.g. with the gravitational 3-body problem, which is known to possess extremely 
complicated bounded orbits. Similarly, the orbits of the harmonic oscillator potential 
V(r) = oo 2 r 2 + 5 are all closed. The content of the celebrated theorem of Bertrand [1] 
is that, under certain technical conditions, the KC and harmonic oscillator potentials 
are the only spherically symmetric ones for which all the bounded trajectories of the 
Hamiltonian iJ(p,q) = p 2 + V(|q|) are periodic. The proof of Bertrand's theorem is 
not difficult, and has actually been included (with various levels of rigor) in several 
textbooks. Nevertheless, this classical result provides a wealth of subtle connections 
with different areas of Physics, and still attracts considerable attention, mainly moti- 
vated by its deep connection with accidental degeneracy and with superintegrability in 
quantum and classical mechanics [2, 3, 4, 5]. 

Some years ago, in a remarkable paper, Perlick [6] showed that Bertrand's theorem 
also arises naturally in General Relativity. Indeed, let us consider a spherically sym- 
metric static spacetime (JA,r)) which is the Lorentzian warped product of the line by 
a Riemannian 3-dimensional (3D) manifold (M,g): 

v = ~dt 2 + g. (1.1) 

Of course, the warping function V G C°°(M) must be strictly positive. It can be easily 
seen that the projection of each timelike geodesic on a constant time leaf M = Mx{t } 
is in fact a trajectory of the natural Hamiltonian H(p, q) = ||p||r* Mo + V(q). This 
relationship between an autonomous natural Hamiltonian flow on a 3D manifold with 
lower bounded potential and the timelike geodesies of a Lorentzian causal manifold 
can be understood as the Lorentzian analog of the introduction of the Jacobi metric 
in classical mechanics [7]. Consequently, by a trajectory in (M.,rj) we shall mean the 
projection of a timelike geodesic onto M . 

Following Perlick, we shall say that a domain B of a smooth (3 + 1)D Lorentzian 
spacetime (A4,rj) (possibly the whole space) is a Bertrand spacetime provided that: 

1. (M.,rj) is spherically symmetric and static, the domain B (possibly minus one or 
two points) is diffeomorphic to a product manifold ]r 1; r 2 [xS 2 x R and the metric 
7] in B takes the form 

77 = dt 2 + g{r) 2 dr 2 + r 2 (dtf 2 + sin 2 $ dip 2 ) . (1.2) 

V(r) 

2. There is a circular (r = constant) trajectory passing through every point of B. 

3. Circular trajectories are stable under small perturbations of the initial conditions. 
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Perlick [6] posed and solved the problem of classifying all Bertrand spacetimes, find- 
ing that the only possibilities are three multiparametric families of spacetimes, called 
hereafter of type I and II±, which are explicitly given by 



Type I 



ds 2 = 



dt 2 dr 2 



+ r 2 (dtf 2 + sin 2 tfd^ 2 ) (1.3) 



Type II ± 



ds 2 = 



G + Vr- 2 + K P 2 (1 + Kr*) 
dt 2 



(1.4) 




-i 



G =F r 2 ( 1 - Dr 2 ± ^/(l - Dr 2 ) 2 - Kr A 



2 1 - Dr 2 ± J (I - Dr 2 ) 2 - Kr 4 




dr 2 + r 2 (dtf 2 + shindy? 2 ) 



where D, G and K are real constants, (3 is a positive rational number, r is a radial 
coordinate restricted to certain interval ]ri,r 2 [, (i?, </?) are spherical coordinates on § 2 
and i G R. 

Bertrand spacetimes are interesting in their own right because of their definition 
in terms of closed orbits. From the point of view of manifold theory, closed geodesies 
have long played a preponderant role in Riemannian geometry [8] . A somewhat similar 
question, that of characterizing all Riemannian manifolds whose geodesies are all closed, 
is still wide open [9, 10]. Periodic trajectories are also central in the theory of dynamical 
systems, where a generic dynamical system with a compact invariant set does have 
many closed orbits [11]. Obstructions do arise, of course, when considering regular 
vector fields whose orbits are all circles. 

The aim of this paper is to prove that, given any (3 + 1)D Bertrand spacetime 
with Lorentzian metric r) (1.1), the term V(r) can always be interpreted as either the 
intrinsic KC or the harmonic oscillator potential on the corresponding 3D Riemannian 
Bertrand space (M, g). Therefore, we do not only have a result on closed trajectories of 
spherically symmetric spaces that extends Bertrand's theorem to 3D curved manifolds, 
but we also recover the key role played by the corresponding KC and harmonic oscillator 
potentials of these 3D spaces. 

In this respect, we stress that Perlick, in a purely relativistic setting, proved that 
KC-like and oscillator-like (3 + 1)D metrics can only be identified provided one consid- 
ers, in this order, Bertrand spacetimes given by the values (3=1 and /3 — 2 [6]. Our 
results show that a proper (non-relativistic) Kepler /oscillator potential can always be 
established for any 3D Bertrand space (M, g), independently of the value of f3, and such 
Kepler /oscillator potentials give all the possible V(r) terms in Perlick's classification 
of Bertrand spacetimes. 

Furthermore, we discuss explicitly some particular Bertrand spaces of interest, such 
as the three classical Riemannian spaces of constant curvature [12] (this case was briefly 
pointed out by Perlick [6]), a Darboux space with non-constant curvature [13] as well 
as the spaces introduced by Iwai and Katayama [14, 15] as a generalization of the 
MIC-Kepler [16, 17, 18] and Taub-NUT metrics [19]. We recall that the Euclidean 
Taub-NUT metric has attracted considerable attention in the physics community (see, 
for instance, [20, 21, 22, 23, 24, 25] and references therein) since, in particular, the 
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relative motion of two monopoles in this metric is asymptotically described by its 
geodesies [26]. Moreover, the reduction of the Euclidean Taub-NUT space by its U(l)- 
symmetry, as performed in the analysis of its geodesies, essentially gives rise to a specific 
Iwai-Katayama space [14] . In this way, all of these different systems are unified within 
our approach in a common Bertrand space framework. 

The article is organized as follows. In section 2 we present the metric structure of 
the 3D Bertrand spaces by considering them within a particular class of spherically 
symmetric Riemannian spaces. In section 3 we identify some physically relevant exam- 
ples of Bertrand spaces. The construction (and interpretation) of the intrinsic KC and 
harmonic oscillator potentials on generic Bertrand spaces is addressed in section 4. In 
section 5, we explicitly study the KC and oscillator potentials for the metrics described 
in section 3. Some comments concerning the generalization of the results here presented 
and their connection with superintegrability properties of the associated Hamiltonian 
flows are drawn in the last section. 



2 Bertrand spaces 

Each of the (Lorentzian) Perlick metrics [6] defines the kinetic energy on the (3 + 1)D 
manifold M.. Alternatively, this construction provides a natural Hamiltonian on the 
3D spherically symmetric space (M, g) with coordinates (r, tp) and with a central 
potential V(r) given by the denominator of the term in dt 2 (1.2). Therefore, according 
to (1.3) and (1.4), hereafter we shall consider the Bertrand spaces defined by 

dr 2 

Type I : ds 2 = ^7^^ + r W + sin 2 d dy, 2 ) (2.1) 

2(l- Dr 2 ± - Dr 2 ) 2 - KrA 
Type II ± : ds 2 = V ^2 ((1 _ ^2)2 _ ^ r 4) " dr " + ^ W + ^ * V)- 

(2.2) 



Consequently, Bertrand spaces are included within the large family of 3D spherically 
symmetric spaces. Their corresponding general metric is usually written in terms of 
a 'conformal' factor depending on a radial coordinate, say f (different from r), and 
sometimes in terms of its logarithm (p = lnf). Hence in order to be able to relate our 
results with others given in the literature we shall consider the following generic metric 
in a 3D Riemannian manifold M written in different coordinate systems: 

ds 2 = /(|q|) 2 dq 2 = f(ff(df 2 + fW) 

= F(p) 2 (dp 2 + dVl 2 ) = g(r) 2 dr 2 + rW (2.3) 

where |q| 2 = q 2 = £? =1 qf and dq 2 = E?=i d % 2 ; /(IqI) = F (p) and #( r ) are 

smooth functions; and df2 2 = d^ 2 + sin 2 ■& dtp 2 is the standard metric on the unit 2D 
sphere S 2 . The relations between the coordinate systems, (q±, q 2 , (fe), (r, 1?, ip), (p, ip) 
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and (r, ip) are given by 

q± — r cos $ 52 = r sin $ cos <p g 3 = f sin i? sin ip 

|q| = f p = lnf F(p)=ff(f) r = F(p) dp = r^g^dr. (2.4) 

The scalar curvature R of the metric (2.3) is always determined by the 'radial' 
(but not necessarily geodesic!) coordinate f <-> p <-> r. This curvature is generically 
non-constant and turns out to be 

A = 2 ,7'(0 2 - 2/(f) (/ ff (f) + 2f- 1 / , (#0) 
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/•» 2 + /'^) 2 -2F(p)F'(^ ,_ 2 ^(r) 2 + 2r^(r)/^(r)-r ^ .^ 



F(p)4 ; V r*g(r) 

where a prime denotes derivative with respect to the corresponding argument, i.e. 
/' = d//df and so on. 

2.1 Spaces of type I 

If we now apply the above results to the Bertrand metrics of type I (2.1) we find that 

9 ^ ~ /VI + Kr 2 dp ~ prVl + Kr 2 ^ 



which yield 



I / ,. x / ,. \ V/3 



p = - In or f = . (2.7) 



Thus the Perlick radial coordinate r is given in terms of the 'conformal' ones by 

2 2 



Hence the metric (2.1) can also be written as 

4 4 



(2.8) 



ds2 = (f-e-Kfey < df2 + fW > = ( e -» -*«*■)' (V + df!2) (2 - 9) 



and its scalar curvature reads 



R = -I {(/3 2 - l)(K 2 f 2 ^ + f" 2 ^) + 2tf(l + 5/3 2 )} 
= -- {(/3 2 - l)(if 2 e 2 ^ + e" 2 ^) + 2tf(l + 5/3 2 )} 
= -4(3/3 2 Kr 2 + /3 2 -l}. (2.10) 
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2.2 Spaces of type II ± 



As far as the Bertrand type II± is concerned, the function g{r) is given by 

_ V2 ( 1 - £>r 2 ± y/(l - Dr 2 ) 2 - Kr 4 \ ^ 

The explicit form of p(r) can be obtained by taking into account that: 

p(r) = f 9 -^dr = fr g -^dr = r f 9 -^dr -/"(/" ^dr) dr 
which yields 

p(r) = rU{r) - J U{r)dr (2.12) 

where 

U(r) = J" ^-dr = Tj^ (l - Dr 2 ± - Dr 2 ) 2 - Kr^ . (2.13) 

The integral f r U(r)dr, although it can be explicitly written, is quite involved, so that 
we omit it and present the results for this family of metrics just in terms of r. 

The scalar curvature turns out to be 

R ^\\h l -^ )+ p 2 (g-gy+i 1 . (2 . 14) 



3 Examples of Bertrand spaces 

In this section we identify within the common framework of the Bertrand spaces (2.1) 
and (2.2) several relevant specific cases that have been introduced in the literature 
by using rather different approaches. In particular, we show that the three classical 
Riemannian spaces of constant curvature [12], a generalization of the Darboux surface 
of type III [13] and the spaces of non-constant curvature studied by Iwai-Katayama [14, 
15] belong to the family of Bertrand spaces. 



3.1 Spaces with constant curvature 

So far we have explicitly shown that Bertrand spaces are generically of non-constant 
curvature. Nevertheless the three classical Riemannian spaces with constant sectional 
curvature k, the spherical (k > 0), Euclidean (k — 0), and hyperbolic spaces (k < 0), 
can be recovered as special instances of Bertrand spaces. The metric of such spaces of 
constant curvature can be collectively written in terms of Poincare coordinates q G M 3 
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(coming from a stereographic projection in R 4 [12]) or in geodesic polar coordinates 
(f, tf, if) as [27, 28] 



sin 2 ( v /fi;f) 
1 + /,•<!-' )-' 1 - ' 7 



d S 2 = ^ ; : dq 2 = df 2 + ^ dfi 2 . (3.1) 



We remark that f is the distance along a minimal geodesic that joins the particle and 
the origin in our Riemannian 3-space; this geodesic distance does not coincide neither 
with f = |q| nor with r (2.4). By taking into account that the conformal function of 
the metric (3.1) is /(f) = 2/(1 + nr 2 ), and by defining 



tanjv^ - ) or r = —= sin (\/ r K f) (3.2) 

we find that the generic metric (2.3) reduces to (3.1). 

This result directly comes from Bertrand spaces (2.1) and (2.2) by setting [6]: 

Type I: (3 = 1 K = -k Type II + : (3 = 2 K = D = k (3.3) 

which yields 



ds 2 = - + r 2 dn 2 (3.4) 



dr 2 

1 — KT 2 

that is, g(r) = (1 — Kr 2 ) -1 / 2 . By applying the change of radial coordinates (3.2), we 
obtain that the latter metric coincides with (3.1). Note also that the scalar curvature 
(2.10) is constant and equal to 6k. 



3.2 Darboux space of type III 

The so-called 2D Darboux spaces were studied by Koenigs [29] and these are the only 
surfaces with non-constant curvature admitting two functionally independent integrals, 
quadratic in the momenta, that commute with the geodesic flow [29, 30, 31, 32, 33]. 
There are four types of such spaces. We shall focus on the Darboux surface of type III, 
T>m, whose metric in terms of isothermal coordinates [34] (u,v) is given by 

ds 2 = e- 2u (l + e u )(du 2 + dv 2 ). (3.5) 
If we define new coordinates (51,52) through 

qi = e ~ u/2 cos(w/2) q 2 = e~ u/2 sm(v/2) (3.6) 
the metric (3.5) is transformed, up to a constant factor, into 

d S 2 = (l + g 2 + g 2 2 )(dg 2 + dg 2 ). (3.7) 

From this expression an iVD generalization of Dm was recently proposed in [13]. Here 
we consider the 3D version given by 

ds 2 — (k + q 2 )dq 2 (3.8) 
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where k is an arbitrary real constant; this metric is clearly of the form (2.3) with 

/(f) = Vfc + f 2 r = F(p) = e p Vk + e 2 P. (3.9) 

Hence 



which gives 



dp Vk 2 + 4r 2 ± k 
g(r) = r — = (3.11) 

V ; dr 2VA; 2 + 4r 2 V ; 



so that the metric (3.8) is transformed into 



2(A; 2 + 4r 2 ) v ; 

In this way we show that Dm is in fact a Bertrand space of type II± (2.2) provided 
that the Perlick parameters are taken as 

(3 = 2 K = D 2 D = -2/k 2 . (3.13) 

The corresponding scalar curvature reads 

2A; + f 2 2k + e 2p k 4 + 2k 2 r 2 - 2r 4 =f k 3 y/k 2 + 4r 2 
R=-6 n -nv, =-6 „ 9 , , =3— ^ • 3.14 

(k + r 2 f (k + e 2 P) 3 r 6 v ' 

3.3 Iwai— Katayama spaces 

These are the spaces underlying the so-called 'multifold Kepler' systems introduced by 
these authors in [15]. Their metric is given by 

ds 2 = r l " 2 (a + br l ' u ) (dr 2 + fW) (3.15) 

where a and b are two real constants, while v is a rational number. This metric is of 
the form (2.3) provided that 

/(f) = r^-\a + br 1 ^) 1 / 2 r = F(p) = (aexp(p/V) + bexp{2p/u)) 1/2 (3.16) 

so that 

f -a±Va 2 + 4br 2 \ 
p = u\n^ j (3.17) 

with b ^ 0. From it we compute the Perlick function g(r): 



dp Va 2 + Abr 2 ± a 

9\ r )= r ^Z.= v / 1 , jf o • ( 3 - 18 ) 

dr Va 2 + 46r 2 



8 



Therefore the Iwai-Katayama metric (3.15) can be written as 

a 2 + Abr 2 

which is proven to belong to the Bertrand family II ± (2.2) under the identifications 

!3 2 = - K = D 2 D = -- . (3.20) 

Hence, the 'multifold Kepler' metric is actually a subcase of the Perlick metric. The 
scalar curvature for these spaces turns out to be 

8ab{u 2 - 1) + Ab 2 [v 2 - I)? 1 '" + a 2 [Av 2 - \)r~ x l v 
~ 2v 2 (a + br 1 l»y 

3a 4 + Qa 2 br 2 + 8b 2 (u 2 - l)r 4 t 3a Va 2 + Abr 2 
Au 2 b 2 r 6 

which acquires a simpler expression once the parameter D (3.20) is plugged in: 



(3.21) 



2(z/ 2 - 1) 1 — Dr 2 =F VI - 2Dr 2 
R ~ v 2 r 2 + 3 • (3 ' 22) 

We stress that the Darboux space of type III with metric (3.12) is a particular 
case of these Iwai-Katayama spaces as it can be recovered by setting b = 1, a = k and 
v — 1/2. Nevertheless a separate study of V m allows us to elaborate on the connections 
with Darboux surfaces and gives rise to a neater analysis of the relationships with KC 
and harmonic oscillator potentials. 



4 KC and oscillator potentials on Bertrand spaces 

It is well known that the KC potential 1/r (r 2 = q 2 = qf + q% + q 2 ) in IR 3 is simply 
a multiple of the (minimal) Green function of the Laplacian A K 3 := J^, whereas 
the harmonic oscillator potential is the inverse of its square. We shall extend this 
prescription, holding on the flat Euclidean space, to a generic curved manifold M with 
metric (2.3). As we shall see, such a prescription is fully consistent with all results 
previously known in the literature. 

The Laplace-Beltrami operator on (M, g) can be written as 

A ^7itM 3g "i (41) 

in terms of an arbitrary set of local coordinates q = (?i,?2,?3)- 111 this formula, as 
usual, g^ denotes the inverse of the matrix tensor ^ in these coordinates and g is the 
determinant of g^. Let us define the radial function C/(|q|), which can be expressed 
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u 



in terms on any of the radial coordinates r, p or r given in (2.4), as the positive 
non-constant solution to the equation 

A M U = on M\{0}. (4.2) 

A short calculation using (2.3) and (2.4) shows that 

1 d / ,dU(f)\ 1 d / , s dU(p)\ 

= -J—iL ( r * iU[r) ) (4.3) 
r 2 g{r) dr \g(r) dr J 

This means that the symmetric Green function U of the Laplace-Beltrami operator in 
M (see [35, 36] and references therein) is given by 

= r_dr_ = /*j£_ = r»M dr (44) 

up to inessential additive and multiplicative constants. Its potential-theoretic interpre- 
tation now leads to define the intrinsic KC potential on the 3D manifold M as 

W KC := a U (4.5) 

where a is an arbitrary constant. The intrinsic harmonic oscillator potential in M is 
defined to be proportional to the inverse square of the KC potential: 

U ho := £ . (4.6) 
In what follows we apply these results to the Perlick metrics (2.1) and (2.2). 

4.1 Type I: intrinsic KC potential 

By considering (2.9) and (4.4) we obtain the Green function U corresponding to the 
Bertrand space of type I with metric (2.1): 



U = (r- p + Kf p ) = (e - ^ + Kef") = -- Vr~ 2 + K. (4.7) 

Then the intrinsic KC and harmonic oscillator potentials on these spaces are defined 
by 

W KC = ~ (r~ P + Kff>) = ~ {e-<» + K&) = ~ Vr^K (4.8) 

U\ 10 = a 7T = a 7T = a — — . (4.9) 

(f-P + KfPf (e-Pp + KePpf r~ 2 + K y J 

Next if we compare these expressions with the warping function appearing in the (3 + 
1)D Bertrand spacetime (1.3), 

V{r) = G + Vr- 2 + K (4.10) 

we find that this is exactly the intrinsic KC potential on the Bertrand space (2.1), 
V = Ukc + G, with a = —(5 and G playing the role of an additive constant of the 
potential. 
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4.2 Type II ± : intrinsic harmonic oscillator potential 



In this case, the Green function U(r) has exactly the expression (2.13). Hence the 
intrinsic KC and harmonic oscillator potentials on the curved spaces of type II± are 
given (only in terms of r) by 

U KC {r) = (l - Dr 2 ± ^J(l- Dr 2 ) 2 - Kr^' 2 (4.11) 

W h o(r) = — a P v (4.12) 

2 1- Dr 2 ± - Dr 2 ) 2 - KrA 

By taking into account the equations (1.4) and (4.12) we conclude that the under- 
lying potential of the (3 + 1)D Bertrand spacetimes of type II±, namely 

V(r) = GTr 2 (l-Dr 2 ± ^(1 - Dr 2 ) 2 - ifr 4 ) ~* (4.13) 

is exactly the intrinsic harmonic oscillator on the corresponding 3D Bertrand space 
(2.2), V = W ho + G, provided that a = =F2//3 2 and G being again an additive constant 
of the potential. 

At this point some comments concerning the above results and those obtained by 
Perlick [6] seem to be pertinent. So far we have obtained all the functions V(r) in the 
families I and II± of (3 + 1)D Bertrand spacetimes as, respectively, the intrinsic KC and 
oscillator potentials on the corresponding 3D Bertrand spaces. This identification holds 
for any value of (3, while in the exhaustive analysis performed by Perlick it is shown 
that relativistic analogues of the KC and oscillator systems can only be obtained by 
considering physical constraints that lead to fixed values (3 = 1 and (3 = 2, respectively, 
of the Bertrand spacetimes (cf. propositions 4 and 5 in [6]). Note that, obviously, both 
results are not in contradiction since the definition of the Kepler /oscillator potentials 
on 3D Bertrand spaces is based on a different (non-relativistic) physical setting that 
does not exclude any value of (3. 



5 Examples of Bertrand spacetimes 

In this section we firstly specialize the intrinsic KC and harmonic oscillator potentials 
for the particular Bertrand spaces described in section 3 and, secondly, we single out 
those Bertrand spacetimes associated to them. 



5.1 Bertrand spacetimes from constant curvature spaces 

Let us consider the classical Riemannian spaces of constant sectional curvature k de- 
scribed in section 3.1 with metric (3.1) or (3.4). The Green function (4.4) is therefore 
given by 

2 f r v 1 
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The definition of the geodesic radial coordinate f (3.2) gives 

tan(v / Kf) r 



2 



In this way the corresponding KC and oscillator potentials are found to be 



(5.2) 



Type I : U KC = -a KT = -ay 'r- 2 - k = -a - — (5.3) 

2r tan{y/ Kr) 

m , 4ar 2 ar 2 tan 2 (-v//«r) . JX 

which are, in this order, particular cases of the Bertrand potential of type I (4.8) and of 
the type II + (4.12) provided that the relations (3.3) have been introduced. Obviously 
we can always add the additive constant G. Notice that, within this framework, we have 
recovered the well known expressions for these potentials (see [27, 28] and references 
therein) . 

Consequently both the KC and harmonic oscillator potentials on these spaces of 
constant curvature lead to particular cases of Bertrand spacetimes; these are 

Type I : d S 2 = - f + + rW (5.5) 

Type II, : ds 2 = — r + + r 2 dQ 2 (5.6) 

+ G -r 2 (2(1 - kt 2 )) l-Kr 2 y 1 

where we have fixed a = — 1 and a = —1/2, respectively. 



5.2 A Bertrand— Darboux spacetime 

In this case the Green function (4.4) of the metric (3.8) or (3.12) turns out to be 



VkTf 2 Vl + ke~ 2 P k±Vk 2 + 4r 2 

U = — — = : = — . (5.7) 

kr k 2kr 

By using this Green function we can define the intrinsic KC and oscillator potentials 
on T>m. However since such a space appears as a particular Bertrand space of type 
II±, only the oscillator potential can be considered in this context; namely 

t tt 7/ OLk 2 r 2 ak 2 2ak 2 r 2 

lypell ± : U ho = ^—~ 2 =— ^ = — — - - (5-8) 

which is reproduced from (4.12) through the identifications (3.13). We recall that this 
was exactly the oscillator potential on T> m worked out in [13] in terms of the radial 
coordinate f = |q|. Thus its corresponding Bertrand spacetime of type II± reads 

, 2 dt 2 k 2 + 2r 2 ± kVk 2 + 4r 2 , 2 2 , - 

ds = l + 0/ , 9 : , 9 x dr 2 + r 2 dVl 2 

G T k 2 r 2 (k 2 + 2r 2 ± k^k 2 + 4r 2 ) 2{k 2 + 4r 2 ) 

(5.9) 

provided that a = =Fl/2. 
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5.3 Bertrand spacetimes from Iwai— Katayama spaces 



The Iwai-Katayama metric (3.15) or (3.19) gives rise to the following Green function 

U = - — \/af- 1 / v + b = - — \TaPpT v + b 
a a 

V2V ( 2 , n , 2 , / o . „ o\ W 



ar 



(a 2 + 2br 2 ± a^a 2 + 4or 2 ) (5.10) 



which allows for the definition of the corresponding intrinsic KC and oscillator po- 
tentials. The results of section 3.3 show that only the oscillator potential can be 
constructed on these spaces in order to obtain a Bertrand spacetime of type II± by 
means of the relations (3.20); explicitly 



Type II ± : W, 



aa 2 aa 2 



ho 



Av 2 {ar- l l v + b) 4u 2 (ae-p/ u + b) 

[a 2 + 2br 2 ± aVa 2 + 4or 2 ) \ (5.11) 



aa 2 r 2 



2v 

These results can be analysed in relation with the so called 'multifold Kepler' po- 
tential Vjk on (3.15) introduced in [15] which is given by 

~2-± 

V lK (f) = r \ U 2 f~ 2 + fi 2 cf-»- 2 + fi 2 drf- 2 ) (5.12) 
a + br» ^ ' 

where /x, c and d are real constants. We remark that this potential is of physical 
relevance as for v — 1 (together with some specific values of the parameters a, b, c and 
d as pointed out in [15]) it comprises the MIC-Kepler problem [16, 17, 18, 37, 38] as 
well as the Taub-NUT one [20, 21, 25, 39, 40, 41]. 

Now notice that the angular momentum of the system L 2 conjugate to df2 2 , ap- 
pearing in the metric (3.15), is a constant of the motion (see e.g. [13]) so that the term 
fi 2 r~ 2 can be reabsorbed within the kinetic energy term since L 2 + /i 2 is obviously an 
integral of motion, /i 2 being an additive constant. Then the potential reduces, up to 
the multiplicative constant /i 2 , to 

V,K(r) = £±£*. (5.13) 
a + br " 

This, in turn, means that VJk can be regarded as a Mobius map on the positive real 
line with variable x = f 1/V = exp(p/V) G M + under the action of the group SL 2 { 

T . dx + c ( d c \ „ r 

x -> Vbc = g • x = - — ■ — where g = G SL 5 

te + a \ o a J 

Next if we add a constant G to the oscillator potential (5.11) we find that 

W ho (r) + G = ^^ (5.14) 

a + or >> 
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provided that 



G = c/a 



ad — be = 



Av 2 ' 



(5.15) 



Note that the second relation is a constraint between the parameters of the potential 
and those of the metric, which arises from the requirement that the Mobius transfor- 
mation be non-degenerate. 

Consequently the 'multifold Kepler' potential Vik (5.13) is actually an intrinsic 
harmonic oscillator (with constant d plus an additive constant c) on the Iwai-Katayama 
spaces studied in section 3.3. 

Finally, we write the resulting metric of the associated Bertrand spacetimes (which 
are of type II±): 



where we have set a = =f2z/ 2 . 

6 Concluding remarks 

In this paper we have explicitly shown that the families of (3 + l)D Bertrand spacetimes 
I and II± can be neatly understood in terms of the 'intrinsic' KC and harmonic oscillator 
potential on the corresponding 3D Bertrand spaces. In this sense, Bertrand spaces 
can be simply described as those manifolds where the classical Bertrand's theorem still 
holds. Moreover, the resulting central potentials cannot be but the KC or the harmonic 
oscillator potentials. Our general description of Bertrand spaces and potentials has 
been also applied to some relevant cases which so far appeared in the literature as 
stemming from different approaches. In particular we have described in a unified way 
the intrinsic oscillator potential on the Darboux III and Iwai-Katayama spaces. As 
a byproduct of these examples we have proven that the so called 'multifold Kepler' 
potential, and therefore the MIC-Kepler and Taub-NUT problems, comes out in a 
natural way within a certain class of (3 + 1)D Bertrand spacetimes of type II±, thus 
suggesting that such potential should be actually regarded as an intrinsic harmonic 
oscillator. 

It is worth recalling that the main algebraic feature of the KC and oscillator po- 
tentials on the iVD Euclidean space is the fact that both Hamiltonian systems are 
maximally superintegrable: this means that both systems admit a maximum number 
(2 iV — 2) of functionally independent and globally defined constants of the motion that 
Poisson-commute with the Hamiltonian. Moreover, it is well known that all bounded 
orbits of any maximally superintegrable Hamiltonian system are periodic. Thus, one 
could conjecture that the KC and oscillator systems on generic Bertrand spaces should 
be maximally superintegrable. In order to prove this statement, the explicit form of 
the 2-3 — 2 = 4 independent integrals of the motion that would commute with these 



ds 2 = 



dt 2 



G T a 2 r 2 (a 2 + 2br 2 ± a^a 2 + Abr 2 ) 




(5.16) 
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Hamiltonian flows should be found. Note that this result is already known in the liter- 
ature for the '//-multifold Kepler systems' with any rational value of v — 1//3 [15]. We 
point out that, again, this problem is quite different to the one consisting in the study 
of the maximal superintegrability of the geodesic motion on the (3 + 1)D Bertrand 
spacetimes which is, in fact, solved (see proposition 6 in [6]) and imposes that (5 — 1, 2. 

We will report on this question elsewhere [42]. For completeness, let us mention that 
this fits naturally into the framework of Hamiltonian systems endowed with an si (2) 
Poisson coalgebra symmetry [43, 44, 45]. It turns out that this coalgebra symmetry 
provides the quasi-maximal superintegrability of any spherically symmetric Hamiltonian 
system defined on a Bertrand spacetime (i.e. three of the above integrals). Therefore, 
in order to achieve the maximal superintegrability of the KC and oscillator potentials 
only one additional integral has to be found by direct methods. The maximal superin- 
tegrability of motion in the particular case of Iwai-Katayama spaces has already been 
established [39, 40, 41], and in fact its (complexified) symmetry algebra is isomorphic 
to that of the KC or harmonic oscillator potentials. Incidentally, it is worth stressing 
that this problem is related to similar questions studied by Koenigs in connection with 
the maximal superintegrability of geodesic flows on surfaces [13, 29, 30, 31, 32, 33]. 

As a final remark, we would like to mention that despite having focused our dis- 
cussion on the (3 + 1)D case, the results here presented can be easily generalized to 
(N + l) dimensions by rewriting the spherical symmetry in terms of an underlying si (2) 
Poisson coalgebra invariance. 
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